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We investigate the spectrum of stochastic gravitational wave background generated by hybrid 
topological defects: domain walls bounded by strings and monopoles connected by strings. Such 
defects typically decay early in the history of the universe, and their mass scale is not subject to the 
constraints imposed by microwave background and millisecond pulsar observations. Nonetheless, the 
intensity of the gravitational wave background from hybrid defects can be quite high at frequencies 
above 10 _8 Hz, and in particular in the frequency range of LIGO, VIRGO and LISA detectors. 
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■ The most interesting cosmological events, from the particle physics point of view, occurred during the first second 
after the big bang. However, the universe became transparent to electromagnetic waves only at the time of decoupling, 
!_i ' tdec ~ 6 x 10 12 (f2/i 2 ) -1 / 2 s, and our ability to observe these events is very limited. In fact, it appears that direct signals 
from sources at t < Is can arrive to us only in the form of gravitational waves. A new generation of gravitational-wave 
detectors will start operating early next century Q , and now is an appropriate time to think about the possible signals 
' we can get from the very early universe. 

A number of potential sources have already been suggested || : quantum fluctuations of the metric during inflation 
HQ , oscillating loops of cosmic string || , and colliding bubbles in a first-order phase transition Q . In this paper we 
£C) . shall investigate the spectrum of gravitational waves generated by transient topological defects, monopoles connected 
by strings and domain walls bounded by strings. 

Among the most significant cosmic events are the phase transitions corresponding to spontaneous symmetry break- 
ings in the early universe. Such transitions are typically accompanied by the formation of topological defects 
Depending on the topology of the symmetry groups involved, the defects can be in the form of surfaces, lines, or points. 
They are called domain walls, strings, and monopoles, respectively. Hybrid defects can be formed in a sequence of 
phase transitions, e.g., the first transition produces monopoles, which get connected by strings at the second phase 
transition. Such defects are transient and eventually decay into relativistic particles. If this happens at a sufficiently 
. early time, the decay products thermalize, and we can see no trace of the defects, except perhaps in the form of 
Q-f gravitational waves ||. In the following, we have adopted units such that c = U = 1. 

Walls bounded by strings. The simplest sequence of phase transitions that results in walls bounded by strings is 
G — » H x Z2 — * H. The first phase transition produces strings, and at the second each string gets attached to a 
C/3 . domain wall. 

The evolution of strings prior to wall formation is identical to that of "ordinary" strings (which do not get attached 
to walls). After a period of overdamped motion, the strings start moving relativistically and approach a scaling 
regime where the characteristic scale of the network is comparable to the horizon. A horizon-size volume at cosmic 
time t contains a few long strings stretching across the volume and a large number of small closed loops. The typical 
distance between long strings and their characteristic curvature radius are both ~ t, but in addition the strings have 
small-scale wiggles of wavelength down to I ~ at, with a <C 1. The typical size of the loops being chopped off the long 
strings is comparable to the scale of the smallest wiggles, I. The loops oscillate and lose their energy by gravitational 
radiation. The lifetime of a loop of size I is r ~ l/TGn, where G is the Newton's constant, n is the mass per unit 
length of string, and r ~ 50 is a numerical coefficient. 

The exact value of the parameter a is not known. Numerical simulations of string evolution [ [To[ suggest that 
a < 10 -3 , while the analysis of gravitational radiation back-reaction indicates that a > TG[i. Thus, we expect a to 
be in the range 10~ 3 > a > TG^i for G^i -C 10~ 5 and a ~ TG^i for G[i > 10~ 5 . We note that for a ~ TG^i the loops 
decay in about one Hubble time from their formation. 

The evolution of the network after domain walls form, having the strings as their boundaries, depends on the 
relative magnitude of the wall tension a, and the force per unit length on a curved string due to the string tension /x, 
F ~ [i/t, where we have assumed that the curvature radius of string is R ~ t. If the walls are formed at t w <C fi/cr, 
then initially they have little effect on string dynamics. The walls become dynamically important at t ~ n/cr, when 
they pull the strings towards one another, and the network breaks into pieces of wall bounded by string. Alternatively, 
if the walls form at t w > /i/tr, then the breakup of the network occurs immediately after the wall formation. We shall 
use the notation t s for the time when the strings start moving relativistically and i* = max{fi/a, t w } for the time of 
the network breakup. The typical size of the pieces is expected to be ~ at*. 
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Gravitational waves emitted by oscillating loops of string during the period t s < t < t* add up to a stochastic 
gravitational- wave background with a nearly flat spectrum fiO-O] , 
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Here, p 9 is the energy density of gravitational waves, p c is the critical density, v is the frequency, 7V"(t) is the effective 
number of spin degrees of freedom in the cosmic plasma at time t, M v is the value of M when waves of frequency v 
were emitted, and f2 r = p r /p c , where p r is the energy density in radiation and relativistic particles. The spectrum 
(Q) extends over a range of frequencies, v min < v < u max , which is determined by the times t s and i». For a rough 
estimate, we can use a simple model which assumes that the loops emit most of their energy in the first few harmonics 
with frequencies v ~ 2/1, where I is the loop's length. To simplify the discussion, we shall assume also that t* is in 
the radiation era and that a ~ TGp (which should be reasonably accurate for superheavy strings with Gp > 10~ 6 ). 
Then, gravitational waves of present frequency v were emitted at 

t„~ 10-^3 (^)"\ (2) 

where /i_ 6 = Gp/10~ e and we have used the approximate form, a oc i 1 / 2 , for the scale factor during the radiation 
era. The cutoff frequencies v m in and v max correspond to the times t* and t s , respectively. In a more realistic model, 
with a full account taken of the radiation spectrum by the loops, the cutoffs would be spread over about two orders 
of magnitude in frequency, while the middle part of the spectrum ([!]) would be unaffected. 

According to Eq.(||), gravitational waves in the frequency range of the LIGO and VIRGO detectors (~ 100 Hz) 
were emitted at i^y ~ 10 _14 /iIgS, while waves in the sensitivity range of the millisecond pulsar (~ 10~ 8 Hz) JlJ] 
were emitted at tup ~ 10 6 pZqS. For ttv < t* < tMP, the gravitational wave background is in the frequency range 
accessible to LIGO/VIRGO, but not to the pulsar. Since in this case the strings decay before decoupling, they leave 
no trace on the microwave radiation, and thus the only constraint on the possible values of Gp comes from the 
nucleosynthesis considerations E3j. 

At the time of nucleosynthesis, the total energy density in gravitational waves, in terms of the critical density, is 
given by 

Kf \ 1/3 / . 
* *nucl \ i / 



0,(nucl) ~ 150G M ^j^j ln^j. (3) 

Here, AC and M nuc i ~ 10 are, respectively, the values of M at t* and at the time of nucleosynthesis, t nuc i ~ 1 s, and 
we have assumed for simplicity that the value of M does not change through the whole period from t s to t* . For the 
standard nucleosynthesis scenario to work, f2 g (nucl) should not exceed 0.054. With AC ~ 100, the resulting bound 
on Gp is 

With the same assumptions, Eq. (jl]) gives 

Q g (u) - 6 x l0- 3 Gph- 2 . (5) 

The logarithm in Eq.(^) can take values from ~ 2 to ~ 100. Consider, for example, strings with Gp ~ 10 -5 and 
walls with symmetry breaking scale r] w ~ 100 GeV. If the symmetry breaking potential for the walls does not have 
very small couplings, then a ~ ryj and t* ~ 1000 s. The time t s when the strings start moving relativistically depends 
on the cosmological scenario. In some models this will be the time when the damping force due to the interaction of 
strings with plasma becomes smaller than the force of tension in the strings, t s ~ (Gp)~ 2 t p ~ 10~ 33 s, where t p is the 
Planck time. Then, ln.(t*/i a ) ~ 80, and the constraint (0) is marginally satisfied. 

The time t s can be greatly increased in inflationary models where the string-forming phase transition occurs during 
inflation, but sufficiently close to its end, so that the strings are not completely inflated away. In this case, t s is 
the time when the typical distance between the strings becomes smaller than the Hubble radius. Strings can also be 
produced in a "pre- heating" transition after inflation Jl6| . The amplitude of scalar field fluctuations at pre-heating can 
reach Planckian values, and superheavy defects can be formed even if the thermalization temperature after inflation 
is very low. The time t s is then shortly after the end of inflation. In some supersymmetric models, certain couplings 
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can be naturally small, and superheavy strings can be formed as late as the electroweak phase transition [jL7j . Then, 
t s ~ t ew ~ 10~ n s. Finally, the effective value of t s can be increased by massive particle annihilations (that is, by a 
decrease of M between t s and t*), resulting in a dilution of gravitational waves from earlier epochs. 

The time is, of course, also model-dependent (through the parameters fi and a), and in general t s and t* do 
not have to be separated by many orders of magnitude. With ln(t»/i s ) ~ 1, the constraint (Q) requires only that 
Gfi < 8 x 10~ 4 , and Eq.(|) gives g (i/) < 5 x 10 6 h 2 . The first version of LIGO interferometer will be able to 
detect a stochastic background with Vl g (v) > 2 x 10~ 6 /i~ 2 after two years of observations, at frequencies ~ 100 Hz and 
with a 90% confidence Q , and thus a detection of gravitational waves from a string- wall network is only marginally 
possible. Advanced LIGO and VIRGO will be sensitive to £l g {v) ^5x 10~ n /i~ 2 , which corresponds to G\i > xlO -8 . 
A similar sensitivity will be achieved by LISA space-based interferometer at v ~ 10~ 3 Hz. 

Monopoles connected by strings. The prototypical sequence of symmetry breakings resulting in monopoles con- 
nected by strings is G — > H x £/(l) — > if. Monopoles formed at the first phase transition get connected by strings 
at the second phase transition. If the monopole-forming phase transition occurs after inflation, then the average 
monopole separation is always smaller than the Hubble radius, and when monopole-antimonopole (MM) pairs get 
connected by strings and start oscillating, they typically dissipate the bulk of their energy to friction in less than a 
Hubble time |§,|l8|. This does not result in any appreciable gravitational- wave background. 

The most interesting scenario, for our purposes, is when monopoles arc formed during inflation (but are not 
completely inflated away). Strings can either be formed later during inflation, or in the post-inflationary epoch. The 
characteristic length scale of strings at formation is then much smaller than the monopole separation; the strings 
connecting monopoles have Brownian shapes, and there is also a large number of closed loops. The evolution of 
strings is initially identical to that of topologically stable strings, without monopoles. At t ~ t s , the strings start 
moving relativistically and generate a gravitational- wave background of intensity (|l]) . In the course of the evolution, 
the characteristic length of strings grows like t and eventually becomes comparable to the monopole separation, so 
that we are left with MM pairs connected by more or less straight strings. We shall call the time when this happens 
t m and the monopole separation at that time, d (note that d ~ t m ). At t > t m , MM pairs oscillate and gradually 
lose their energy by gravitational radiation. 

In the above discussion we have assumed implicitely that t m > t s . In this case, using the inequalities d ~ t m > 
t s > (Gfi)~ 2 t p , it is easily verified that 

lid S> m, (6) 

where m is the monopole mass. Then the motion of monopoles is relativistic, with a typical Lorentz factor 

7 ~ fj,d/m > 1. (7) 

If strings are formed during inflation, soon after the monopoles, then the strings connecting MM pairs can be nearly 
straight. In this case, the strings do not go through a period of relativistic evolution, and no gravitational radiation 
is produced prior to t m . But unless these defects are formed very close to the end of inflation, the inequality (fij) will 
still be satisfied. We shall assume this to be the case. 

Oscillating MM pairs lose their energy by gravitational radiation and by friction due to scattering of plasma 
particles off the monopoles Q. To estimate the rate of gravitational energy loss, we calculated the power and 
the spectrum of gravitational waves radiated by an oscillating pair connected by a straight string. Details of this 
calculation will be given elsewhere [ po[ , while here we shall only summarize the relevant features of the radiation 
spectrum. Most of the energy of an oscillating MM pair is radiated in the frequency range 

cT 1 < < 7 2 eT 1 (8) 

with a spectrum 

d£ g /dtdv ps AG^v- 1 . (9) 
At higher frequencies, d£ g /dtdv oc v~ 2 . The total radiation power is pl| 

£ g =fGfi 2 , (10) 

where 

fps81n7 (11) 
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is a numerical coefficient taking values in the range 10 < f < 10 3 , depending on the Lorentz factor 7. 

A monopole-antimonopole pair connected by a straight string is, of course, a very special configuration, and one 
could be concerned that the radiation spectrum in this case may be very different from that for a generic configuration. 
To address this issue, we calculated the spectrum of electromagnetic radiation emitted by a pair of oscillating, equal 
and opposite charges connected by a straight string. We found that the resulting spectrum is qualitatively similar 
to that in the generic case plj . This suggests that the same may be true in the case of gravitational radiation, but 
of course a detailed analysis is required before any definite conclusion can be drawn. Here, we shall assume that the 
spectrum (|^) is representative of the generic case [ p3[ . 

The frictional energy loss is of the order 

tf ~ T 2 , (12) 
where T is the temperature, and it is easily verified that £f £ g for t m 3> t s . Hence, the lifetime of a pair is 

t - fid/£ g ~ d/TGfj,. (13) 

The number density of AIM pairs at t ~ t m ~ d is njvf(t m ) ~ d~ 3 , and at later times n^i(t) ~ (id) -3 / 2 (t m < 
t < r) . The spectral power of the gravitational wave background emitted by the oscillating pairs per Hubble time at 
time t is 

Q.f{v) ~ 4G[i 2 tn M (t)/ Pc (t) ~ 120(Gfi) 2 (t/d)^ 2 , (14) 

where p c (t) ~ l/30Gi 2 is the critical density. The highest power comes from the radiation emitted at t ~ t. The 
frequency range (||) at that time corresponds to the present range 

v* < v < 7 2 z/*, (15) 

where 

v n ~ d-\r/t eq ) l / 2 z- q \ (16) 

i eg is the time of equal matter and radiation densities, and z eq is the corresponding redshift. The spectral power in 
the range (|l5| ) is [compare with Eq.((|)] 

Sl B (v) ~ 120f- 3 / 2 (G^) 1/2 (AA/AA r ) 1/3 ^r- (17) 

For i/<i/», the main contribution to the spectrum is made at the time ~ (f/i/*) 2 r, when the frequency 
corresponding to v at the present time is ~ d~~ 1 . The resulting spectral power is given by Eq.(|l7|) with an additional 
factor (yjv^. This part of the spectrum extends to the minimum frequency v min ~ (rG/i) 1 / 2 ^*. 

Radiation emitted prior to t nuc i ~ Is should satisfy the nucleosynthesis constraint. Assuming that r < t nuc i, the 
total energy density in gravitational waves at t nuc i is 

n g (nucl) ~ 30(G/x/f ) 1/2 (A/; uc i/A/;) 1/3 < 0.05, (18) 

where we have used Eq. (|ll]) . Combining this with Jl7|), we have 

0,(zy) < lO- 5 ?- 1 ^ 2 (19) 



For the frequency range ( |l5| ) to overlap with the sensitivity range of LIGO/ VIRGO, it is necessary that d > 
2 x 10 _13 (rG^) _1 cm, which implies T < 400 (assuming that m < m p , where m p is the Planck mass), and thus 
^g( v ) ^2x 10 _8 /i -2 . This is below the sensitivity of the first version of LIGO but may be within the reach of the 
advanced detectors. 

If the strings did have a period of rclativistic evolution, then the spectrum includes another flat region with spectral 
power (H) at v > (fG^)" 1 / 2 ^,. This power is greater than ( [l7|) if G\i > f -3 which, for f ~ 10 3 , includes all 
interesting values of G/z. The analysis of this part of the spectrum is identical to that for walls bounded by strings, 
with i» replaced by t m . As before, the advanced detectors will have sufficient sensitivity in the case of strings with 
G/x < xl0~ 8 , provided that the flat portion of the spectrum includes v ~ 100Hz. 
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